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DOCKET NO- 3351-064 

DIRECTION FINDING METHOD AND SYSTEM USING DIGITAL 
DIRECTIONAL CORRELATORS 

Field of the Invention 

[0011 The present invention relates generally to direction finding systems, and in particular, 
to a digital direction finding method and system utilizing a revolving antenna array. 

Background of the Invention 

[002] The direction finding of electromagnetic transmission signals has been practiced for 
many years and has been achieved using a wide variety of antenna configurations and signal 
processing techniques. One such antenna configuration is a linear array of receiving antennas. 
Each receiving antenna has its own receiver circuit, and they are configured to simultaneously 
receive the transmission signal. At long distances, the radiating source appears as a point source, 
and the transmission signal from the point source appears as a planar wave-fi-ont to the linear 
array of a relatively small length. 

[003] The simultaneously received signals firom each receiving antenna are processed to 
determine the phase difference of the received signals between a pair of the receiving antennas 
(baseline). The measured angle of arrival of the transmission signal is a function of the distances 
between the receiving antennas in the baseline, the frequency of the transmission signal, and the 
measured phase difference. Thus, if the distance between the receiving antennas of the baseline 
is known, and the frequency and the phase differences are obtained through measurements, then 
the angle of arrival of the transmission signal relative to the array can be determined. If the array 
position in inertial space is known, then the angular position of the emitter can be determined in 
any coordinate system. 

[004] A system of this type is called an interferometer direction finder. In general, 
interferometer direction finders are capable of high precision angle of arrival measurements, but 
the angle measurement has multiple ambiguities caused by the periodic nature of the phase 
difference. 
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[005] Conventionally, multiple interferometer baselines must be used to resolve the 
ambiguities, known as multiple baseline interferometry. In multiple baseline interferometry, 
different observations of the angle of arrival by several baselines are correlated to remove the 
ambiguities of the angle of arrival. Disadvantages of this approach are that it requires the space 
to install multiple antennas, preferably in a coplanar geometry, and it requires extra antenna 
elements and receiver electronics to process the multiple baseline signals. Moreover, to obtain 
precision angle measurements for VHF (very high) and UHF (ultra high) frequencies (50 MHz to 
200 MHz), a long interferometer baseline is needed. The greater the end-to-end antenna spacing, 
the more ambiguities are introduced. The greater the number of ambiguities, the greater the 
number of intermediate antennas w^hich must be added. This represents a large hardware cost 
and is undesirable, especially, if the antenna array is to be installed in a limited space. 

[006] The direction finding problem is also difficult in practical application because the 
fi-equency of the incoming signal typically is not known before its reception at the direction 
finding system. Li fact, the direction finding system may be required to operate over a very 
broad frequency band and an incoming signal may lie anywhere within the band. An approach 
in such a case is to implement the direction finding system as a channelized device and employ 
digital channelizers to divide the broad frequency band at each antenna into a number of plural 
fi-equency channels (each narrower than the broad frequency band) to be processed individually. 
The main problem with this approach is that a relatively large number of digital channelizers are 
required, i.e., one for each antenna. 

[007] The direction finding system can be combined with a parameter encoder system 
which measures various parameters including frequency, pulse width, etc. of the transmission 
signal. Alternatively, the signals received at the receiving antennas can be converted through a 
Fourier transformafion, followed by a phase-sequence analysis to determine the transmission 
signal frequency. Again, additional extemal devices or processing circuits are necessary. 

[008] Moreover, if the direction finding antennas are broad-band, the angular coverage of 
the arrays is typically 90*=* (spatially) and four such arrays are needed to cover fiiUy 360"". 

[009] Thus, there is a need for direction finding method and system capable of determining 
the angle of arrival with high accuracy using as few antennas as possible , and requiring no 
dedicated circuit for the signal frequency determination. 
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Summary of the Invention 

[010] It is an object of the present invention to provide direction finding method and system 
which are capable of determining the angle of arrival of a wave-front with high accuracy using a 
few antenna elements, and requiring no dedicated circuit for the signal frequency determination. 
The method and system of the present invention are especially suitable for detecting distant 
radiating sources in the RF frequency bands. 

[Oil] It is a further object of the present invention to provide direction finding method and 
system which determines both frequency and direction of arrival by virtue of the cross- 
correlation processing technique. 

[012] These and other objects of the present invention are achieved by a method of 
detecting a wave-front, comprising the steps of measuring a first relational phase between signals 
induced by the wave-front in first and second antenna elements of an antenna array at a first 
spatial position. The antenna array is then moved to a second spatial position and a second 
relational phase is measured again. This process is repeated until all ambiguities have been 
resolved and the angle of arrival and frequency of the wave-front can be determined. 

[013] The above and other objects of the present invention are further achieved by a 
direction finding method using a non-stationary antenna array. In accordance with the method, a 
set of antenna array spatial positions and relational phases between signals induced in the 
antenna array is predicted for each potential value of at least one parameter of a wave-front. 
Then, while the antenna array is moved through a plurality of spatial positions, actual relational 
phase at each of the spatial positions is measured to generate a set of measured spatial positions 
and relational phases. Finally, the measured set is compared with each of the predicted sets to 
find a best-matclung predicted set and thus the ambiguities are resolved. 

[014] The above direction finding method of the invention can be implemented by 
hardware, software or both. In a hardware implementation, a directional correlation system of 
the invention comprises an antenna array; a steering mechanism for moving the antenna array; a 
receiver for measuring a relational phase therebetween antenna elements of the antenna array; 
and a correlator for determining a parameter of a wave-front incident to the antenna array. 

[015] The correlator includes a memory for storing predicted relational phases, that are 
expected to be found in the anteima array, in a plurality of memory cells each being associated 
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with a potential value of the parameter to be determined. An encoder is provided to shift 
contents of the memory cells corresponding to the movement of the antenna array. A set of 
comparators are coupled to receive outputs of the receiver and access selected memory cells for 
matching the measured relational phase fed by the receiver with the content of the selected 
memory cells. A set of counters are provided to count a number of matches for each memory 
cell during the movement of the antenna array. A threshold detector is coupled to outputs of the 
counters for determining the memory cell the content of which matches best with the measured 
relational phase, thereby determining the parameter of the wave- front. 

[016] In accordance with an aspect of the invention, the antenna array may comprise more 
than one baseline. 

[017] In accordance with another aspect of the invention, the steering mechanism rotates 
the antenna array in either planar or conical surfaces. 

[018] In the former case, the system can determine the frequency and direction of arrival of 
the plane wave-front whose normal lies in the plane of rotation of the antenna arrays. If the 
normal to the wave-front is inclined to the plane of rotation, then the system can determine the 
angle of arrival of the projection of the normal on the plane of rotation, and one of the frequency 
of the wave- front and the angle of inclination, provided that the other is known. 

[019] In the later case, frequency, azimuth and elevation angles can all be uniquely 
determined. 

[020] In accordance with a further aspect of the invention, the correlator is a multi- 
dimensional correlator. In one instance, the multi-dimensional correlator comprises a column (or 
row) of uni-dimensional correlators. In another instance, the multi-dimensional correlator 
comprises a matrix of uni-dimensional correlators. 

[021] Still other objects and advantages of the present invention will become readily 
apparent to those skilled in the art from the following detailed description, wherein the preferred 
embodiments of the invention are shown and described, simply by way of illustration of the best 
mode contemplated of carrying out the invention. As will be realized, the invention is capable of 
other and different embodiments, and its several details are capable of modifications in various 
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obvious respects, all without departing from the invention. Accordingly, the drawings and 
description thereof are to be regarded as illustrative in nature, and not as restrictive. 

Brief Description of the Drawings 

[022] The present invention is illustrated by way of example, and not by limitation, in the 
figures of the accompanying drawings, wherein elements having the same reference numeral 
designations represent like elements throughout and wherein: 

[023] Figure 1 is basic geometry showing the orientation of an electromagnetic wave-front 
relative to an antenna baseline; 

[024] Figure 2 is a schematic diagram of an exemplary interferometric phase measurement 
circuit which can be used in the system of the present invention. 

[025] Figure 3 is rotating array geometry showing the orientation of the electromagnetic 
wave-front of Figure 1 relative to the antenna baseline when the antenna baseline rotates for an 
angular interval; 

[026] Figure 4 is a conceptual processor for use in the system of the present invention; 

[027] Figure 5 is an implementation of the conceptual processor shown in Figure 4; 

[028] Figure 6 is a schematic diagram of an embodiment of the system of the present 
invention using the implementation illustrated in Figure 5; 

[029] Figures 7-10 are graphs showing the auto-correlation fiinction for various numbers of 
quantization bits; 

[030] Figure 1 1 is a schematic diagram of another embodiment of the system of the present 
invention having two revolving baselines; 

[031] Figure 12 is a schematic diagram of a digital cross-correlator for use with the 
embodiment shown in Figure 1 1 ; 

[032] Figure 13 is a schematic diagram of still another embodiment of the system of the 
present invention using a two-dimensional cross-correlator; 
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[033] Figure 14 is three-dimensionally rotating array geometry showing two baselines of 
the system of the invention in three dimensional movement; 

[034] Figure 15 is a schematic diagram of yet another embodiment of the system of the 
present invention using a three-dimensional cross-correlator; and 

[035] Figure 16 is a high-level block diagram illustrating an exemplary computer system 
with which the present invention can be implemented. 

Best Mode for Carrying Out the Invention 

[036] Direction finding method and system according to the present invention are 
described. In the following detailed description, for purposes of explanation, numerous specific 
details are set forth in order to provide a thorough understanding of the present invention. It will 
be apparent, however, that the present invention may be practiced without these specific details. 
In other instances, well-known structures and devices are shown in block diagram form in order 
to simplify the drawing. 

1. Mathematical Background 

[037] Figure 1 is basic geometry showing the orientation of an electromagnetic wave-front 
101 having a normal 100. Antenna array 102 comprises two receiving antennas 104, 106 spaced 
apart a known distance d to form a baseline of antenna array 102. The relative electrical phase, 
or phase difference, y/, of the voltages induced in receiving antennas 104, 106 by the wave-front 
101 is given by the following equation: 

2nd ^ . ^ 
y/- cos^Jsm^ (1) 

where the relative electrical phase in radians, 

d is the end-to-end spacing between receiving antennas 104, 106, 
X is the wavelength of wave-front 101, 

<!> is the elevation angle, and 6 is the azimuth angle of normal 100 with respect to the 
azimuth plane which is also the horizontal plane in this particular example. 
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[038] For sake of simplicity, the antenna array 102 is shown in Figure 1 to have a boresight 
direction 108 which is coincident with axis x. The boresight direction of an antenna array is the 
direction of strongest radiation in the polar pattern. As depicted in Figure 1 , azimuth angle 0 is 
measured between the boresight direction 108 and the radiation direction (projection 100) in the 
horizontal plane. Likewise, elevation angle ^ is measured between the normal 100 and its 
projection 109 in the horizontal plane. 

[039] Equation (1) can be rewritten in a normalized form to obtain equation (2). 

y/ = 27ts{[ + A:)cos ^ sin 0 (2) 

where d = sAq , and 

— = (l + /:)— or / = (l + /c:)/o,/is the frequency of wave-front 101. 

[040J Thus "5" represents the end-to-end spacing of receiving antennas 104, 106 in 
wavelengths of a reference wavelength Ao, and (1 + k) represents the frequency of 
electromagnetic wave-front 101 in units of the reference frequency fo which corresponds to wave 
length ylo. 

[041] If phase difference if/ is digitized (quantized) to obtain an integer N satisfying 
0 < AT < 2" -1 , where w is an integer representing the number of quantization bits, then equation 
(2) can be rewritten to obtain equation (3) 

N = Int[(2"s{l + fc)cos^ sin0)Modulo2" J (3) 

where A'^ is the integer part of the term in brackets. 

[042] Since sin^can have both positive and negative values, and is defined as a positive 
integer, negative remainders in the modulo operation must be increased by the addition of 2" -1 . 
Equations (4.1) and (4.2) describe the value of for Pranging from zero to 360 spatial degrees, 
and ^ ranging from -90 degrees to +90 degrees. 

N = Int[[2"s{l-hk)cos^sm0)Modulo2"l forsine>0 (4.1) 
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;V = 2" - 1 - Int[(2" s{l + k)cos <z>|sin 0\)Modulo2" ], for sinG < 0 (4.2) 

[043] The foregoing discussion is true with the conventional stationary interferometer. The 
present invention, in contrast, involves revolving antenna arrays to resolve the ambiguities 
encountered in determination of the frequency and direction of arrival of a wave-front. Figure 3 
is rotating array geometry showing the orientation of electromagnetic wave- front 101 relative to 
antenna array 102 which is at an angular interval a with respect to the zero azimuth direction. 
Then, azimuth angle 0 in equations (1) through (4) must be replaced by {0-a). Equation (4) 
can be written in fimctional form to obtain equation (5) 

N = g[(l + ^)cos ^ s\n{e - a)] (5) 

[044] Equation (5) is easily modifiable for other types of planar movement, such as linear 
movement or combination of linear and rotating movements. It should also be understood that 
the azimuth plane is not necessarily a horizontal plane (parallel to ground surface). 

[045] Equation (5) states that the number (or phase difference i//) is related to the values of 
k (or frequency J), elevation angle ^, azimuth angle 0 and angular interval a (or the spatial 
position of antenna array 102). Equation (5) is ambiguous in the sense that there is more than 
one set of values of k, ^, 0 and a which correspond to each value of N. 

2. Rotating the Interferometer to Determine the Wave- front Parameters 

[046] The primary object of a direction finding system is to deduce the wave-front 
parameters (0, ^, k) from measurements on the system. Let U = (1 + A:)cos^ then equation (5) 
becomes 

N = g[Usmi0-a)] (6) 

[047] The discussion will be restricted in the following section to a single distant pulsed 
radiating source which produces wave-front 101 at the interferometer or antenna array 102. The 
discussion will then be extended to a multiplicity of radiating sources. 
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• • • . 

[048] As antenna array 102 rotates, angular interval a and digitized phase difference are 
measured. If there were no ambiguities, two measurements would be enough to determine 
azimuth angle ^and intermediate variable U (i.e. two equations and two unknowns). However, 
because each pair of values {N, d) may correspond to many pairs of values of {6, U), a great 
many pairs of values (AT, a) must be collected to determine uniquely actual values (0, U), From 
equation (6), the set of values \N,a\ can be predicted for each (^^,{7_^). Thus, there will be a 

particular set 



N, a, 



N, a. 



corresponding to each {6^,Uy)y where r = 1, 2,..., n. Each Nr will range from 0 to 2" - 1, and 
each Qv will range from 0° to 360^. 

[049] A processor 400 (shown in Figure 4) can be conceived for which each {p^^Uy) 
(where X = 1, 2..., /; and y = 1, 2..., J) a set [A^,Q^];c,;; is predicted from equation (6) and stored. As 
antenna array 102 rotates, the direction finding system measures the electrical phase difference (// 
and outputs a quantized value N corresponding to each antenna spatial position, a, for which a 
signal from the radiating source has been received. The measured set of values [A^,arl„^^, is 
compared against each of the stored set of values [N.a^^^^^j . The particular set of stored values 
[A^,Qr]^^ for which the best match is obtained will correspond to the wave-front parameters 
{Ojc^Uy). A best match is defined as that set [N,a\f^^^ which has the most members which 
match the set of measured values [A'^jCir]^^^, . 
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[050] Since f/ = (1 + k)cos^ , k or frequency / must be known to determine elevation angle 
^ and vice-versa. Thus, the described processor 400 can determine azimuth angle 6 and 
frequency ^, or azimuth angle 0 and elevation angle ^. The interferometer may be tuned to a 
specific frequency band such that the system will respond only to that frequency band. Then 
elevation angle ^ and azimuth angle ^may be determined. Alternatively, antenna array 102 may 
be so designed and so oriented such that the system will respond to signals arriving at a 
particular elevation angle In this case, frequency ^ and azimuth angle ^can be determined. 

3. Svstem Design 

[OSl] One way of implementing processor 400 is through the use of a shift register structure 
500 illustrated in Figure 5. Consider the case where a multiplicity of pulsed emitters are all at 
the same (or nearly the same) elevation angle and the same (or nearly the same) frequency. That 
is, the emitters in this particular example all have the same (or nearly the same) value of U. It is 
desired to determine their azimuthal direction through the use of the rotating interferometer. 

[052] Equation (6) shows that the set [N,a] can be predicted for an emitter at any azimuth 
angle ^ since is established. Register structure 500 of Figure 5 shows how this information 
may be stored in processor 400. In Figure 5, each cell of register 500 represents 0.5"^ in azimuth. 
With the interferometer boresight direction at 0° as shown, the value of N for an emitter at any 
azimuth angle is stored in register 500 within the limits of resolution determined by the quantum 
steps of 0.5*^. Thus the transition from one value of N to the next value of TV can be represented 
to within ± 0.25'' for register 500 of Figure 5. 

[053] As the interferometer rotates, the code will be shifted from cell to cell in concert with 
the rotation. Every time the boresight direction of antenna array 102 changes by 0.5°, each value 
stored in a cell will be shifted into the adjacent cell in the direction of the boresight change. If 
the value contained in a particular cell is accessed as shown in Figure 5, then the set of values 
produced by the shifting action described previously, coupled with the antenna boresight 
position, will describe the set [N,a] for the particular register cell (^a) that is accessed. 
Therefore, if a multiplicity of register cells are accessed, a multiplicity of sets [N,a] can be 
generated, each corresponding to a different azimuth angle ft 
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[054] Figure 6 shows a possible embodiment of the system of the present invention using 
register structure 500 shown in Figure 5. 

[055] Phase difference y/ of the rotating interferometer 600 is measured by receiver 602. 
Receiver 602 can be any device available in the art which is capable of measuring the phase 
difference between signals induced in receiving antennas 104, 106. An exemplary principle 
schema of receiver 602 is illustrated in Fig. 2. Receiver 602 may further comprise a digitizer (not 
shown) which produces a digital code A^r or giO^) of the measured phase difference, wherein 0^ 
represents the azimuth angular position of the emitter. 

[056] Antenna array 102 of rotating interferometer 600 is rotated by a steering mechanism 
604 which can be any device available in the art. Rotating interferometer 600 may also include 
an encoder 610 for shifting the content of shift register 500 according to the rotation of antenna 
array 102. Preferably, encoder 610 is coupled with steering mechanism 604. 

[057] Shift register 500 has the entire code pattern stored in it for 360° of azimuth. Shift 
register 500 is accessed at various cells 502 corresponding to specific azimuth angular positions 
of the emitter (e.g. ft, ft-i,...)- The received phase difference or code g{B^ is compared against 
the stored codes g{Oi) for all accessed cells 502. For example, if cells 502 were accessed every 
2^, there would be 180 lines coming from register 500 feeding 180 AND gates 606. The code 
for each received pulse would be compared simultaneously with the code at each AND gate 606. 
Wherever a match occurs between the received code g{9^ and the stored codes g(^), the 
corresponding counter 608 is incremented by one count. Apparently, the counters may be 
configured as reverse counters. If at least one pulse were received from an emitter for each of 
the 0.5° shift register cells over 360°, then the counter at the output of the AND gate 
corresponding to the emitter's azimuth angular position would have a total count of 720. Li this 
particular example, AND gates 606 fimction as comparators. 

[058] It should be understood that once a match occurs at a particular AND gate 606, that 
AND gate must be locked until the antenna array position changes, and data in shift register 500 
also changes by at least one cell shift. This is important because counters 608 is simply 
recording the fact that for a given orientation of the antenna array, the received code matched the 
stored (reference) code. Repeated pulses from the emitter when the antenna array position has 
not changed merely affirm this condition, add no new information and must be locked out. For 
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this purpose, a blocking circuits 612, which can be any device and connected to each AND gates 
606 in any manner available in the art, is provided. For example, blocking circuit 612 may feed 
output signals to an AND gate 606 to optionally enable/disable the AND gate 606 as desired. 

[059] Once one (or more) of counters 608 reaches a specified level (at, or near the 
maximum count, i.e., about 75%), a threshold detector 614 is actuated and an emitter is declared 
present at the azimuth angular position corresponding to the counter 608 in question. Because 
the phase difference code is ambiguous, more than one AND gate 606 will have a match between 
the received code from a particular emitter pulse and the stored codes. Thus some counters 608 
will have counted to some value less than the count in the counter corresponding to the emitter's 
azimuth angular position. 

[060] It is important to note that the term "shift registers" and "AND gates" used in the 
foregoing description are in general meant to be generic terms, and are accurate descriptors only 
for the case where « = 1 (1 bit). For the case where w > 1 , both terms, "shift registers" and 
"AND gates" mean multibit devices, which perform the indicated operations on multibit words. 

1061] The process described above is a cross-correlation between the received code and the 
codes stored in shift register 500. It should be understood that shift register 500 is just one 
hardware implementation of processor 400 and other arrangements of processor 400 are also 
available. For example, any of shift register 500, AND gates 606, counters 608, blocking circuit 
612 and threshold detector 614 can be implemented as a piece of software executed in a 
compliant computer envirorunent. In a pure software embodiment, processor 400 may be 
implemented in form of a multi-dimensional table or database which stores for each potential 
value of azimuth angle B a set \Nr,c\tored' The measured set of values [A^,<3r]^^^^ is compared 
against each of the stored set of values [A'^,^],^^^^^ to find the best match stored set. The table or 
database can be accommodated within an ordinary computer system such as an exemplary 
computer system 160 shown in Figure 16. 

[062] Computer system 160 includes a bus 162 or other communication mechanism for 
communicating information, and a processing unit 164 coupled with bus 162 for processing 
information. Computer system 160 also includes a main memory 166, such as a random access 
memory (RAM) or other dynamic storage device, coupled to bus 162 for storing information and 
instructions to be executed by processing unit 164. Main memory 166 may also be used for 
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storing temporary variables or other intermediate information during execution of instructions to 
be executed by processing unit 164. Computer system 160 further includes a read only memory 
(ROM) 168 or other static storage device coupled to bus 162 for storing static information and 
instructions for processing imit 164. A storage device 170, such as a magnetic disk or optical 
disk, is provided and coupled to bus 162 for storing information and instructions. 

[063] Computer system 160 may be coupled via bus 162 to a display 172, such as a cathode 
ray tube (CRT) or a flat panel display, for displaying information to a computer user. An input 
device 174, including alphanumeric and other keys, is coupled to bus 162 for communicating 
information and command selections to processing unit 164. Another type of user input device is 
cursor control 176, such as a mouse, a trackball, or cursor direction keys for communicating 
direction information and command selections to processing unit 164 and for controlling cursor 
movement on display 172. This input device typically has two degrees of freedom in two axes, a 
first axis and a second axis allowing the device to specify positions in a plane. 

[064] According to one embodiment of the invention, a variety of information and services 
are provided by computer system 160 in response to processing unit 164 executing sequences of 
instructions contained in main memory 166. Such instructions may be read into main memory 
166 from another computer-readable medium, such as storage device 170. However, the 
computer-readable medium is not limited to devices such as storage device 170. For example, 
the computer-readable medium may include a floppy disk, a flexible disk, hard disk, magnetic 
tape, or any other magnetic medium, a CD-ROM, DVD-ROM, or any other optical medium, 
punch cards, paper tape, or any other physical medium with patterns of holes, a RAM, a PROM, 
an EPROM, a FLASH-EPROM, or any other memory chip or cartridge, a carrier wave embodied 
in an electrical, electromagnetic, infrared, or optical signal, or any other medium from which a 
computer can read. Execution of the sequences of instructions contained in main memory 166 
causes processing unit 164 to perform the process steps described below. In altemative 
embodiments, hard-wired circuitry may be used in place of or in combination with computer 
software instructions to implement the invention. Thus, embodiments of the inven 

[065] Computer system 160 also includes a communication interface 178 coupled to bus 
162. Communication interface 178 provides a two-way data communication as is known. For 
example, communication interface 178 may be an integrated services digital network (ISDN) 
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card, a digital subscriber line (DSL) card, a cable modem, or a telephone modem to provide a 
data communication connection to a corresponding type of communication line. As another 
example, communication interface 1 78 may be a local area network (LAN) card to provide a data 
communication connection tion are not limited to any specific combination of hardware circuitry 
and software, to a compatible LAN. Wireless links may also be implemented. In any such 
implementation, communication interface 1 78 sends and receives, through a data communication 
connection 180, electrical, electromagnetic or optical signals which carry digital data streams 
representing various types of information. Of a particular note, the communications through 
interface 178 may permit transmission or receipt. 

[066] In an embodiment of the invention, interface 178 of computer system 160 is coupled 
to received measured data set [A^,a^l„eflj from receiver 602 and steering mechanism 604. Then, a 
dedicated software can be executed to match the received code with the stored codes. Thus, the 
fields of the table or database fimction as the cells or memory cells of register 500. 

[067] In another embodiment, measured data set [^,cir]„g^^ observed by receiver 602 and 
steering mechanism 604 may be temporarily stored in a removable storage media, such as a 
floppy disk of CDROM, and transferred later to computer system 160 for processing. 

[068] One important factor in a cross-correlation process is the main lobe-to-side-lobe ratio. 
[069] Figure 8 shows the auto-correlation function for n — 2, and Figure 9 shows the auto- 
correlation function for « = 4 bits. In each case the side lobe level is on the order of — of the 

main lobe level. Since the code is not a perfectly random function, and seems to repeat itself 
every 360 electrical degrees (repetition is not exact because of the factor sin(^-a) in equation 
5), the near-in side lobes exceed the average value by a goodly margin, as can be seen in e.g. 
Figure 7. 

[070] Figure 10 shows the behavior of the main lobe and the first few side lobes as a 
function of n. The main lobe becomes narrower as n increases. 

[071] If more than one emitter were present, then the counter corresponding to each 
emitter's azimuth angular location would register a maximum (or near maximum) count. The 
average count in other counters 608 (those not corresponding to each emitter spatial angle) 
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would be — where M is the number of emitters. In order to differentiate the main lobes from 

the total average side lobe level, the main lobe should be on the order of 4 times the average side 
lobe level. Therefore 

2" 

— < 4 or AT < 2""' 
M 

Then a system with az = 5 would have a finding capacity of about 8 emitters. 

4. Lowering the Auto-correlation side lobes 

[0721 An upper bound on n for the interferometer is established by the signal-to-noise 
(S/N) ratio, and receiver measurement errors. The number of emitters that can be handled 
before the system saturates is in effect limited by the S/N ratio which results in reduced 
sensitivity. In order to lower the auto-correlation side lobes, a second antenna baseline and a 
second receiver (RF processor) are added as shown in Figure 11. 

[073] In Figure 1 1, a rotating interferometer system or antenna array 1 10 is formed with 
two baselines 114 and 115 having end-to~end spacing d] and d2, respectively. Although 
baselines 114 and 115 may be formed from four individual antenna elements, it is sufficient 
to use only three antenna elements 111, 112 and 113. Signals from the shared antenna 
element 1 1 1 are supplied to both RF processors 1 16 and 1 17 corresponding to baselines 1 14 
and 115, respectively. RF processors 116 and 117 are of a type similar to receiver 602 
shown in Fig. 6, and need not be described in detail. 

[074] In one embodiment, both baselines 114 and 115 are driven by a single steering 
mechanism 604 (not shown) and a single encoder 610 is needed to feed the angular positions 
of antenna array 110 to the subsequent correlator circuit. In another embodiment (not 
shown), baselines 114 and 115 are driven separately by two steering mechanisms similar to 
steering mechanism 604, and two encoders similar to encoder 610 are needed to feed the 
angular positions of antenna array 1 10 to the subsequent correlator circuit. 

[075] The subsequent correlator circuit or digital cross-correlator 120 is schematically 
depicted in Figure 12. Digital cross-correlator 120 includes two channels 121 and 122. Each 
channel is almost identical to the correlator circuit shown in Figure 6. That is, each channel 
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comprises a shift register (123, 124) and a plurality of comparators (125, 126). Chamiels 
121, 120 share a common set of AND gates 128, counters 608 and threshold detector 614. 

[0761 Shift registers 123, 124 store the entire code pattern for 360° of azimuth for 
baselines 114, 115, respectively. Each of shift register 123, 124 is accessed at various cells 
corresponding to specific azimuth angular positions of the emitter or emitters. The received 
phase difference received fi-om RF processors 116, 117 are compared, by comparators 125, 
126, against the stored codes of the corresponding shift register 123, 124 for all accessed 
cells. A match must be obtained on each channel before a count is registered in a given 
counter 608. This is ensured by supplying the comparison results of each of comparators 125 
and its counterpart 126 to one of AND gates 128. The counter 608 corresponding to the 
AND gate 128 is incremented by one count only when a match occurs in each of the 
comparator 125 and its coxmterpart 126. Threshold detector 614 outputs the actual azimuth 
angular position of the emitter in the manner described with respect to Figure 6. In essence, 
correlation is performed in two charmels in parallel. 

[077] Assuming each channel 121, 122 has the same number of quantization «. Then, 
the total number of bits available becomes In and the number of emitters that ran be handled 
(based on previous reasoning) before the system saturates is given by 

;V = 2'"-' =(2"-')' (7) 

By adding a third antenna or a second baseline, the capacity of the system is significantly 
increased. For example, a system with n = 5 would now have a finding capacity of about 256 
emitters. That is, 256 emitters can be simultaneously detected by rotating interferometer 
system 110 with adequate accuracy. 

[078] It has also been determined through computer modeling that the two baselines 
114, 115 should not be equal in length, nor should they be co-linear to obtain the lowest side- 
lobe-to-main-lobe ratios. The ratio d\ld2 is preferably in the range between 1 and 2, and the 
angle p is preferably in the range between 40 and 70 degrees. The concept can be extended 
to more than two baselines by adding more antennas and more receivers, and performing 
correlation in parallel as described for the two baseline case. 
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5. Determination of both t/and 9 

[079] Up to this point the discussion has been restricted to finding the direction of 
emitters at a particular frequency / and elevation angle if> (the same value of CO- To 
determine both U and 6, a two dimensional correlation must be performed, e.g. in a system 
130 shown in Figure 13. Each of the cross-correlators 120 shown in Figure 13 is identical to 
that shown in Figure 12. However shift registers 123, 124 in each cross-correlator 120 have 
different codes stored in them corresponding to the value of U for the particular correlator 
with which the shift register pair is associated. Shift registers 123, 124 of each cross- 
correlator 120 are connected to their respective RF processors 116, 117 via a corresponding 
input bus 131, 132. The outputs of threshold detectors 614 of cross-correlator 120 are 
consolidated into an output bus 135. For each emitter, only the counter in the cross- 
correlator corresponding to the actual 6 and XJ values, respectively, of the emitter will have a 
maximum count. 

[080] Again, it should be understood that more than two baselines can be used in the 
above two dimensional interferometer system by adding more antennas, receivers and shift 
registers to each of the cross-correlators. Likewise, the two dimensional interferometer 
system shown in Figure 13 does not necessarily have two or more baselines. One base line 
should also be sufficient to simultaneously determine both 6 and U values of one or more 
emitters. As discussed above, adding more baselines and correlation channels, however, 
helps to enhance the system accuracy and finding capability. 

[081] Again, the two dimensional interferometer system can be realized in a software 
implementation in the manner described above. 

6. Determination of Elevation 6 and Frequency f 

[082] It has been shown in Sections 1 and 2 that if the interferometer is rotated such that 
the baseline is always in a plane (i.e. azimuth plane), the phase difference is given by 
equation (8): 

= Imi). + /:)cos <l>^ %\n{e^ - a) (8) 
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where ^ and 9^ are the emitter elevation and azimuth angles, respectively. It was also noted 
in Section 3 that either frequency / or elevation angle ^> can be obtained if the other is known, 
but both cannot be obtained along with emitter azimuth angle ^by rotating the interferometer 
in a plane. 

[083] If the interferometer is rotated such that its baseline (or baselines) lies on a conical 
surface (or surfaces) as shown in Figure 14, then the relative phase difference between the 
two antenna elements of the baseline is given by equation (9). 

i// = 2;zs(l + A:)[cos cos ip^ cos{0^ ~ ) + sin sin ] (9) 

where is the elevation angle of the baseline and 0a is the azimuth angle of the baseline 
with respect to the azimuth plane. Apparently, steering mechanism or mechanisms 604 in 
this case must be able to drive the antenna array in three dimensional movement. 

[084] The subscripts 1, and 2 in Figure 14 refer to the two interferometers which are 
used jointly in the same manner as the system of Figure 11. In general the spacing between 
the two antenna elements of each interferometer is different (as in Figure 11). The phase 
difference code for each interferometer is given by (10.1) and (10.2) 

=Int^''{\■hk)s^[cos^^^cos^^cos{0^-0„^)'^sm^^^sm^ (10.1) 
If /«^{ }, <0,then N, =2" -l + /nr{ }, 
A^2 =/«/{2"(l + A:)yJcos<*„2COS^^cos(^^ -(9^J+sin^^2sin<*^]M^ (10.2) 

If /«/{ }^ <0,then =2"-l + //i/{ }^ 

[085] Azimuth and elevation angles and frequency would be determined by cross- 
correlators similar to that shown in Figure 12. The system could be organized as shown in 
Figure 15. In Figure 13, the correlator circuit includes a column 134 of a plurality of cross- 
correlator of the type shown in Figure 12. In Figure 15, the correlator circuit includes a 
plurality of columns 134 each corresponding to a specific value of one of the parameters of 
the emitters, i.e. elevation angle ^. Each cross-correlator 120 in a column 134 corresponds to 
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a specific value of one of the remaining parameters of the emitters, i.e. fi-equency / (or k). 
Thus, each cross-correlator 120 has a code pattern stored in its shift registers corresponding to 
the particular fi-equency and elevation angle [k^.^j] for which the cross-correlator is 

matched. The code is obtained firom equations (10.1) emd (10.2) for each shift register 
respectively. Thus, three dimensional correlation is performed in this case. 

[086] If frequency /is known.(e.g. by virtue of using super-heterodyne receivers), then 
each colunm 134 will have only one cross-correlator 120, and only a single row 139 is 
required. If elevation angle ^ is known (e.g. by using directional antenna elements responsive 
only to radiation incident to the antenna array at a particular elevation angle), then only a 
single column 134 is required. In each instance, system 150 is simplified to have the form of 
system 130. 

[087] Again, more than two baselines can be used in the above three dimensional 
interferometer system by adding more antennas, receivers and shift registers to each of the 
cross-correlators. Likewise, one base line should be sufficient to simultaneously determine 
all three parameters 6 , ^ and / of one or more emitters. 

[088] Again, the three dimensional interferometer system can be realized in a software 
implementation in the manner described above, 

[089] It should be understood that the system and method of the invention is not limited 
to determination of azimuth, elevation and firequency of a wave-front. Instead, the system of 
the invention is capable of determining any parameter of a radiating source which can be 
detected by antenna elements of a revolving antenna array as long as the parameter can be 
expressed as a function of spatial position of the antenna array and phase difference induced 
in the antenna elements. 

[090] Moreover, azimuth, elevation and frequency play interchangeable roles, 
especially in the software implementation of the invention. For example, the multi- 
dimensional table or database which stores a vast number of [A^,a] sets each corresponding 
to a particular set of (0, J) values may be indexed in any desired order, e.g. first by 
elevation angle ^, then by azimuth angle 0, and finally by frequency / Alternatively, the 
multi-dimensional table or database may be indexed to present a vast number of ^,/| sets. 
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each corresponding to a particular pair of (iV, a). A comparison software module then 
compare each stored {N, a) pair with the data received from the system's receivers and 
encoders to highlight the [0, (p, f] sets corresponding to the matched (A^, a) pair. The 
particular [6, ^, /| set which is most frequently highlighted will be outputted as actual 
parameters of the wave-front. Apparently, the parameters 0^ ^, and / are treated equally in 
this embodiment. 

[091] The method and system of the present invention have many advantages over the 
known direction finding technique. For example, a single antenna array, instead of four 
arrays, is sufficient to cover 360^ spatially. Moreover, it is sufficient to rotate the system's 
antenna array through an angular sector on the order of 60 to 90 degrees. Upon completion 
of the rotation, frequency and direction of arrival are available on all emitters within the 
angular sector. The information is also available continuously as the antenna array rotates. 
The system finding capacity, i.e. the number of emitters for which the system can measure 
frequency and direction of arrival, can be predetermined based on the precision of the 
electrical phase difference measurement, and the number of antenna elements utilized. 

[092] In addition, the tectmique of the present invention involves far less antenna and 
receiver hardware than the fixed (non-rotating) interferometer system. Though, the system of 
the invention cannot determine direction of arrival on a single pulse, after a single rotation of 
the antenna array, the frequency and direction (both azimuth and elevation) on all emitters 
within the system's operating band and sensitivity will be available. The complexity of the 
system, in the hardware embodiments, lies in the digital circuitry which performs the 
necessary cross-correlation for direction of arrival and frequency determination. The greater 
the system resolution in frequency and direction of arrival, the greater the digital complexity 
(in terms of number of digital circuits). In contrast, the conventional interferometer must add 
antenna and receiver channels to achieve better angular resolution. 

[093] While there have been described and illustrated specific embodiments of the 
invention, it will be clear that variations in the details of the embodiments specifically 
illustrated and described may be made without departing from the true spirit and scope of the 
invention as defined in the appended claims. 
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